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Abstract: We study the collider implication of a neutral resonance which decays to several 
diboson final states such as W + W~, ZZ, and Z 7 via a minimal set of effective operators. 
We consider both CP-even and CP-odd bosonic states with spin 0, 1, or 2. The produc¬ 
tion cross sections for the bosonic resonance states are obtained with the effective operators 
involving gluons (and quarks), and the branching fractions are obtained with the operators 
responsible for the interactions with electroweak gauge bosons. We demonstrate that each 
scenario allows for a broad parameter space which could accommodate the recently-reported 
intriguing excesses in the ATLAS diboson final states, and discuss how the CP states and 
spin information of the resonance can be extracted at the LHC run II. 


Keywords: LHC, diboson, resonance, effective operators 



Contents 


1 Introduction 1 

2 Spin-0 resonances 2 

3 Spin-1 resonances 6 

4 Spin-2 resonances 12 

5 Kinematic Correlations in the Diboson Final State 15 

6 Summary 17 

A Decay widths 18 

A.l CP-even scalar 18 

A.2 CP-odd scalar 19 

A.3 CP-even tensor 19 


1 Introduction 

Recently, the ATLAS collaboration has reported some excesses in searches for diboson reso¬ 
nances in the highly boosted final states with W + W~, W±Z and ZZ at the 8 TeV LHC with 
20.3 fb -1 [1]. They have adapted boosted techniques to tag hadronically decaying W and Z 
gauge bosons, which strongly suppress the QCD dijet backgrounds. All three excesses emerge 
at around 2 TeV in the invariant mass distribution formed by two W- or Z-tagged fat jets. 
The CMS collaboration also sees a moderate excess at the similar location in all hadronic 
channel [2, 3]. In response to the tantalizing experimental observations, several papers have 
already appeared taking this phenomenon as a new physics signature [4-30]. 

A typical recipe for a new physics model to explain the above-mentioned excesses is the 
introduction of two new heavy states: a charged particle and a neutral particle. The former 
takes care of the W±Z channel while the latter does the other two channels. However, given 
the fact that a large fraction of events belong to all three channels, it may be a reasonable 
attempt to fit the data only with a single new heavy resonance. As a matter of fact, Allanach, 
Gripaios and Sutherland recently investigated the diboson resonances in this direction: they 
basically introduced a likelihood function for the true signal in the W + W~, W±Z, and ZZ 
channels and found that the maximum likelihood has zero events in the W±Z channel [15]. If 
this observation were true, the ATLAS data would indicate a single neutral bosonic resonat¬ 
ing particle rather than two, which show up in all three channels due to detector effects and 
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misidentification of W Y and Z bosons. We also note that in the single particle interpretation 
coincidence of the resonances at 2 TeV in the three channels can be naturally understood. 
Keeping this minimality and simplicity of the single particle interpretation, we further inves¬ 
tigate the possible classification of neutral resonances by considering different spins and CP 
states in an effective field theory approach including a set of operators for each case. 

Our philosophy is basically the bottom-up approach, invoking a minimal set of effective 
operators that may be responsible for the W + W~ and ZZ signals. As no spin informa¬ 
tion is available, we extensively consider spin-0, spin- 1 , and spin -2 resonances. Symmetries 
of the relevant operators also induce potential signals in different final states, encouraging 
experimental collaborations to look into the related channels for consistency. 

In the following three sections, we examine scalar, vector, and tensor resonances in turn, 
focusing on viable parameter scans in conjunction with production cross sections and partial 
decay widths of the resonance at hand. In Section 5, we briefly make comments on kinematic 
correlations among the final state particles to extract spin, CP states, and coupling informa¬ 
tion of the resonance of interest and the proposed interactions. Section 6 is reserved for a 
summary. 

2 Spin-0 resonances 

In our new physics interpretation, the resonance particle decays into two bosons so that the 
resonance itself should be a bosonic state with an integer spin. In this section, we begin 
with considering a spin-0 resonance and study the effects of its CP states with corresponding 
effective operators. 

A CP-even scalar resonance (henceforth denoted as S) in diboson channel could be well- 
parameterized by the following interaction Lagrangian: 

C s = ~\s(s x F^F y ^ + s 2 F™F w ^ + s :i G%G a ^+ £ s/m///), ( 2 . 1 ) 

/ 

where F Y V and F^ Y denote the field strength tensors for usual U(l)y and SU(2)w gauge 
bosons before the electroweak symmetry breaking while denotes the SU(3) C gluon field 
strength tensor with the color index a = 1, 2, • • • 8. 1 The strengths of the above couplings are 
parametrized by si, s 2 , and s 3 , respectively for gauge bosons and s/ for fermions. 

A tiny flavor non-diagonal interaction would lead un-acceptable flavor changing neutral 
current (FCNC) effects so that we naturally expect that the coefficients s/ are negligibly 
small or strictly flavor diagonal. The first generation quarks, u and d, could have the largest 
contribution to the production of the scalar at the LHC but they are suppressed by a small 
factor ~ nif / A. Furthermore, the coefficient s/ can be forbidden by a global symmetry when 
the singlet scalar is promoted to a complex scalar field T with Re(T) = S. The couplings to 

X A scalar particle such as gravi-scalar or radion [31-33] potentially provides diboson resonance and may 
have other signatures [34-36]. However, we found that the width of 2 TeV radion is unacceptably big to 
account for the ATLAS anomaly. 
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the gauge bosons in the form of Eq. (2.1), however, are still obtained due to SM anomalies. 
A similar argument can be applied to the CP-odd spin-0 resonance. Considering all these, we 
would take the gluon fusion as the dominant production mechanism for the scalar resonance 
and neglect the production by diquark. 

Without loss of generality, we take S 3 = 1 by redefining A. The other coefficients, s 1 
and s 2 , for U(l)y and SU(2)w gauge kinetic terms, are redefined as relative strengths to s 3 . 
From the interactions in Eq.(2.1), we obtain the partial decay widths of S into 77 , Z 7 , ZZ, 
W + W~, and gg as 


S-y 7 | 2 r7lg 

47tA 2 


r s (77) = 

P s(ZZ) = yjl - 4xf (1 - 4x| + 6(x§) 2 ) , 

r ' S {zi) = (1 - x s z )\ _ 

Y S (W + W~) = '^V l- 4 ^ (1 - + 6(40 2 ) 


[T s (gg) = ^ 


?l 2m s 


s 77 = si cos 2 9w + S 2 sin 2 9\y , 
Szz = S 2 cos 2 9w + Si sill 2 9\y , 
sz'y = (s 2 - si)sin2(9jy , (2.2) 

swtv = 2 s2 

s 99 = S 3 , 


where ms and 9\y denote the mass of CP-even scalar S and the Weinberg angle. Here and 
henceforth, we define the mass squared ratio of a heavy SM boson i (Z, W, or h ) to a 
resonance R as 




TO, 


TO 


R 


(2.3) 


Obviously, in this parametrization, S is produced via gluon fusion followed by the decays 
into the above final states. Of potential experimental constraints, the two following conditions 
should be settled to be in the right “ball park” with respect to the recent ATLAS data: 


• the total decay width should be within ~ 10 % of the mass of the resonance [ 1 ] , 

• the signal production cross section should be as sizable as order of several fb [15]. 

In general, the single production cross section of a narrow resonance is proportional to the 
total decay width of the decaying particle. Therefore, demanding a sizable production cross 
section with a (relatively) smaller total decay width is not a trivial task. We remark that 
as discussed in the literature, reported excesses in all three diboson final states {W + W~, 
W±Z and ZZ) are not independent of one another, and the data in one channel may be 
contaminated by the data in the other channels due to detector effects. As we mentioned in 
introduction, in Ref. [15], authors performed a general analysis of new physics interpretations 
of the recent ATLAS diboson excesses by computing a likelihood function for the true signal 
in the W + W~, W±Z, and ZZ channels. They found that the maximum likelihood has zero 
events in the W±Z channel, i.e., one could fit the data in all three channels with a single 
neutral resonance in the final state with W + W~ and ZZ. The likelihood is sufficiently flat 
and the required cross section (for 95% C.L.) is in the following range [15]: 

0(4 - 8 ) fb < a ■ BR(W + W~) + a ■ BR(ZZ) < 0(20 - 24) fb , (2.4) 
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Figure 1: The left panel shows production cross sections (in fb) of the CP-even scalar resonance in 
the final states with W + W~ + ZZ (red solid curves) and Z 7 (blue dashed curves). The black dotted 
curves represent the constraint of the total decay width, T/mg. The dark yellow-shaded region is 
excluded by the dijet search and the light green-shaded region represents the allowed parameter space. 
The right panel shows branching fractions of the spin-0 resonance (Z'y, ZZ, W + W~, and gg by red 
dotted, blue dot-dashed, magenta dashed, and black solid curves, respectively) as a function of si. 


where a is the single production cross section of the resonance. For our analysis with the case 
of the CP-even scalar, we first fix the mass of the scalar resonance, ms, to 2 TeV, and then 
compute the signal cross section, cr(pp —> S —> W + W~ + ZZ), by varying three parameters, 
A, si and S 2 (s 3 =l). We find that in the majority of parameter space, the consistency (gauge 
invariance and Lorentz invariance) of the model predicts a large branching fraction into the 
diphoton final state. In particular, when two parameters have the same sign (i.e., S 1 S 2 > 0), 
diphoton rate (oc |s 77 | 2 = |sicos 2 0 w + S 2 sin 2 0 iy| 2 ) turns out to be too large so that the 
model is severely constrained by current data at the 8 TeV LHC [37, 38] . 

Interestingly enough, the opposite case with S 1 S 2 < 0 constraint provides a way to reduce 
the diphoton rate as clear from Eq. (2.2). Especially, the condition of si ~ — tan 2 0 p/S 2 
gives s 77 ~ 0 thus vanishingly small diphoton final state, and the condition fixes the relative 
branching fractions as follows: 

BR(W + W~) : BR(ZZ) : BR(Z 7 ) : BR(gg) 

]si | 2 _ I Si 1 2 cos 2 26V . ]si | 2 . 1 , 2 

4 tan 4 8 sin 4 6\y 4 tan 2 6\y 

We take this relation for illustration and calculate the relevant cross sections in the two 
dimensional parameter space of A vs. si, although other relations can be straightforwardly 
analyzed. For the relevant data analysis (and remaining analyses throughout this paper), 
we employ Monte Carlo event generators CalcHEP [39] and MadGraph5_aMC@NL0 [40]. In 
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Fig. 1, we show production cross sections (in fb) of the CP-even scalar resonance in the 
final states with W + W~ + ZZ (red solid curves) and Z 7 (blue dashed curves). Contours of 
T/ms are shown by black-dotted curves. The corresponding branching fractions are shown 
in the right panel as a function of si. The dijet resonance searches provide constraints (at 
95% C.L.) on the parameters, which are shown by the dark yellow-shaded region [41, 42]. 
Combining all constraints together, the allowed parameter space represented by the light 
green-shaded region might accommodate the diboson excesses. We remark that the exact 
relation of si = — tan 2 9\ys 2 is not required, and any minor deviation from this relation would 
be easily allowed as long as the associated diphoton rate is below the current limit [37, 38]. 

Speaking of CP-odd spin-0 case, a pseudo-scalar or axion-like scalar (denoted as A) can 
couple to the SM gauge bosons through anomalies. The gauge interactions are parametrized 
in a way similar to the CP-even scalar case with one of the field strength tensors replaced by 
a dual field strength tensor: 

L a = -jA( a] F^F Y ^ + a 2 F™F Wpv + a 3 G^G pv ) , (2.6) 


where the dual field strength tensors are defined as, for example, F^ v = \e^upaF 5 pcr , and 
prefactors 01 , a 2 , and a 3 denote the coupling constants which can be determined by anomalies 
for a global symmetry. For instance, a*/A = CiOti/{Sn/ a) (i = 1,2,3) with Ja being the 
breaking scale of a global U( 1) and q = K) Q 9 o^Gi(^a) where q a is the global U( 1) charge of 
a heavy fermion and tGi{ r a ) is the Dynkin index for a representation r a under the SM gauge 
group Gi [43, 44], 

The total decay width of the pseudo-scalar resonance [43] is given by the sum of partial 
decay widths into 77 , Z 7 , ZZ, W + W ~, and gg: 


'r A hi) = ^\c^, 

r ^( z b) = £ydcz 7 | 2 (l -®z) » 

< Va{ZZ) = ^\c zz \ 2 (l-Axif 2 , 
r a{W+W~) = ^\cww\ 2 (l ~ 4x^-) 3/2 , 
jA(gg) = 


c 77 = aq cos 2 9\y + sin 2 9w , 

cz 7 = (a 2 - ai) sin (29 w ), 

czz = a 2 cos 2 9 W + ai sin 2 9 W , (2.7) 

cww = 2 a 2 , 

c 99 = a 3 , 


where xf is defined in Eq. (2.3). The case with the CP-odd scalar has similarities compared 
to the case with the CP-even scalar in the sense that the corresponding branching fractions 
are similar along with associated coefficients, and also we require aia 2 < 0 to suppress the 
diphoton rate. Like the CP-even scalar case we simply choose a 2 = — ai/tan 2 f?vv, and 
demonstrate the resulting parameter scans in Fig. 2. We observe that all contours are similar 
to those in Fig. 1, except for the scale of A due to a larger cross section for the CP-odd scalar. 

A couple of comments should be made here. Speaking of the unitarity bound for scalar 
resonances first, we observe that in both CP-even and CP-odd cases, the spin-0 resonance 
couples only to transverse modes of SM gauge bosons. Then, the unitarity cutoff can be 
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Figure 2: Same as in Fig. 1 but for CP-odd scalar resonance. Branching fractions are similar to 
those for the CP-even scalar, replacing s, by at- 



just read from the coefficients of the effective operators, namely, of order max(A/sj) and 
max (A/a,;) in CP-even and -odd cases, respectively, by power counting. Thus, as shown in 
Fig. 1 and Fig. 2, the unitarity cutoff is > A ~ 10 TeV, which is consistent with the effective 
interactions with a TeV-scale resonance. Second, we find that in both CP-even and CP-odd 
cases, the Z 7 production cross section is about 1-3 fb at the 8 TeV in the allowed parameter 
space. The current experimental data tells that the (95% C.L.) upper bound on the Z 7 
production in the dilepton channel is given up to the resonance mass of 1.6 TeV while the 
higher mass reach is limited by statistics [45] . Nevertheless, we expect that the corresponding 
limit for the 2 TeV resonance would be comparable to the result at the resonance mass of 
1.6 TeV in Z 7 searches or below the existing limit. Therefore, a(pp —>■ Z 7 ) = 0(1 ) fb is still 
allowed for the 2 TeV and this channel would rather provide an interesting consistency check 
for ZZ and W + W~ excesses. As shown in the right panel of Fig. 1, BR(Z^f) is comparable 
to BR(ZZ ), and one cannot turn off BR(Z 7 ), as it would also eliminate the signal. In other 
words, if diboson excesses turned out to be the real signal with a CP-even or -odd scalar, 
observation of an excess in the Z 7 channel would corroborate the case. 

3 Spin-1 resonances 

For spin -1 resonances, we consider an extra U(l)x gauge symmetry that is realized by the 
Stueckelberg mechanism. Then, the would-be Goldstone boson ax ensures the gauge invari¬ 
ance of the effective action. 

First, imposing the SM gauge symmetry and U( l)x, we have the dimension-4 interaction 
Lagrangian between the U(l)x gauge boson and the quarks and/or gauge bosons in the SM 
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given as follows: 


C-DA = -gxZ^idPL + c R P R )q - \eF x u F x ^ - (ir,D»a x (H^D p H) + c.c.) , (3.1) 

where the covariant derivative is defined as D p a x = d lt a X — (Jx with a x being the Stueck- 
elberg axion, g x is the Z l gauge coupling, and c R = cl (c r = —cl ) for CP-even (CP-odd) Z'. 
These dimension-4 interactions correspond to diquark couplings, gauge kinetic mixing and 
mass mixing in order. 

We keep the dimension-4 diquark coupling to Z l in Eq. (3.1) as a production mechanism, 
while the lepton couplings are suppressed as in leptophobic Z' models [46] . The gauge kinetic 
mixing with e / 0 leads to T{Z' Zh) = Y{Z' —> WW) due to the SM gauge symmetry. As 
the Zh channel is strongly constrained by the LHC bound, a(pp — > Z') x BR(Z / — > Zh) < 7fb 
[47], which is significantly lower than the required value, ~ 10 fb, for explaining the ALTAS 
diboson excesses [27]. Moreover, no ZZ decay is induced from the dimension-4 operators. 
Therefore, we do not consider the possibility of a sizable mass mixing with Z' any more taking 
e<l, Instead, we consider novel effective interactions for Z' containing the ZZ decay mode, 
coming from dimension-6 operators, as will be discussed below. The last term with real r\ in 
Eq. (3.1), which is a CP-even interaction, should be highly suppressed, because of potential 
Z l decays into Z R h or WlW r . If 77 is purely imaginary, namely, the last term in Eq. (3.1) 
is equivalent to a CP-odd operator (d^D^a^W H up to a total derivative, thus leading to 
(d^Z'^H^ H , but a vanishing on-shell decay amplitude squared for Z' —> H. 

It is noteworthy that the interactions of a vector isospin triplet W' to the SM electroweak 
gauge bosons can be introduced by a similar dimension-4 operator in the effective theory such 
as W' p cr a H [15], which mixes the extra gauge boson with the SM massive gauge bosons. 
In this case, the ATLAS diboson excesses can be explained by the W±Z channel, provided 
that the charged spin-1 resonance is produced via quark annihilation at the LHC [15]. In 
our work, we do not investigate the potential of the charged resonance as mentioned earlier 
because the dibosonic decays of a neutral resonance suffice to explain the ATLAS diboson 
excesses within current experimental errors. 

Moving onto higher dimensional operators, we enumerate CP-even dimension-6 operators 
as follows [48]: 

Crm = Y 2 D^a x [i(D^H^F x H + c.c.] + g D^a x [(D u H)^F X H + c.c.] 

+g D‘*a x [i(D t 'H)iF% / H + c.c.] + ^ D^a x [(D 1 ' il) t F%H + c.c.] 

d^D^ax {biFj^'P* + b 2 F™F w P" + b 3 G pa G , (3.2) 

where A is of order the mass of extra heavy fermions, and a* (z = 1, 2, 3,4) and bi (i = 1, 2,3) 
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parametrize the coupling strengths. The CP-odd counterparts of dimension-6 interactions are 


tm = jr 2 D ll a x {i(D I 'H)^F^H + c.c.] + g D ll a x [(D u H)^F^ U H + c.c.] 
+^D‘ i a x [i(D’'H)*F%,H + c.c.] + ^D»a x [(D''H)*F$H + c.c.] 

+Y d^D^a x (j)\Fj rT F Yp<J + b 2 F™F wpa + b 3 G p(T G pt7 ), (3.3) 

where a* ( i = 1, 2,3,4) and bi ( i = 1, 2, 3) parametrize the coupling strengths. We comment on 
the dimension -6 operators composed of one held strength tensor for Z l and two held strength 
tensors for the SM gauge bosons: Tr( F p x F Xv F UIL ) for CP-even operators and Tr (F YX F Xl/ F Ufl ) 
for CP-odd operators with F pu = F Y V . F^ Y , G pu . First of all, the CP-odd operators can be 
rewritten as F p u F uX F Xp = F pu F Xu F pX , which is the same as F^F Xv F" x = -F*F Xu F pX , and 
as a result, we get F p v F vX F xp = 0. Likewise, the CP-even operators can be also rewritten 
as F p v F vX F x>l = F pi/ F Xu F pX . Then, using the identity of F Xp F uX = — F a pF a P, we get 
F Yv F uX F XlJi = —^F^ p F a pF a P = 0. Therefore, the dimension -6 operators composed of gauge 
held strength tensors only are identically zero so that we do not consider them in our analysis. 

Given the above observations, the Z l gauge boson decays only by symmetry breaking 
terms given in Cd& or £d 6 - When it comes to the production modes for the spin-1 resonance, 
we henceforth assume that it is produced by diquark couplings and ignore the gauge kinetic 
mixing and mass mixing. The effective cubic interactions for Z' coming from Cd 6 and Cdq 
are obtained as shown below: 

£cp—even = ^ (a 1 m z Z' / Z'^F* + a 2 d v h Z^F^j 

~y(\ axmze^Z^(d p W"a ~ daW^ p ) + a A SFhZ’^W* - d u W^ 

+ m X? Z l{ ~ \a 3 e^W~(d p W+ - d a W+) + ia A W-"(d,W+ - d u W+) + c.c.) 

+Y &tZ l { hF P<T pYp(T + b2F^F wp ° + b 3 G pa G p ^ , (3.4) 

£ C P-odd = y ijM rnzZ^Z'^F^ + a 2 d v h Z ,p F Y ^ 

-Y^mzZ^Z'^ 9 ^ - d v W Y ) + X -a^ pa d v hZ\^d p Wl - d CT W p 3 )) 

+ ^pz; t ( - a 3 W-(& i W +v - 8 U W +P ) + l -i~a A e^W-{d p W Y - 8 a W+) + c.c.) 

+ Y 9>XZ l (biF Y pa F Yp ° + b 2 F%F wp ° + hG p<J G p °) , (3.5) 


where the U(l) x gauge coupling is absorbed into ai and ai, and so on. 



After the electroweak symmetry breaking (EWSB) and dropping the terms with the 
divergence of Z 7 , the effective CP-even interactions for Z 7 are 


£ C P-even = «i e pUp ° Z' p Z u F pa + Z' p Z v (d p Z a - d a Z p ) 

+(K 2 e^ pa z' p w-{d p w+ - d a w+) + ik 2 z' p w- u (d p w+ - d„w+) + c.c.) 

+J- Z' p d v hF pv + ^ Z'v&'h {d p z v - d v z p ) (3.6) 

where F pv is the photon field strength tensor and k 2 , k 2 are related to other parameters by 
gauge invariance as 

k 2 = ——— (k\ sin 0w 4- ki cos 6\y ), k 2 = — —(^-(^3 sinf?vv + k^cosOw)- (3.7) 
mz k / A 


We note that the effective triple gauge interactions with Z l in the above effective Lagrangian 
are the generalized Chern-Simons terms that are generated by extra heavy fermions [44, 48]. 
Using the effective action above, we obtain the partial decay rates of the spin-1 resonance [44] 
into Zy, ZZ, W + W~, /iy, hZ, and qq, respectively as: 


r z'(zz) = i^(i-44') ; Z 


1 + x z 


^ 



1 — x 


r*<»z)-fiE&(i- 

X (2 + sf'(xf - x z zf + 2xf(xf + 3x§') - (4sf + 3®§')) 
k rz'(99) = Sa ^-(l + 2xf )(l -4xf) 1 , 


(3.8) 


where is defined in Eq. (2.3). On the other hand, the effective CP-odd interactions for 
Z 7 become 


£cP-odd = ot\Z' p Z v F pu + diZ^Z^Z, - d v z p ) 

+ [a 2 Z' p W~^d p W+ - d v W+) + ia 2 e pvpa Z' ll W-{d p W+ - d ff W+) + c.c) 

+ y £ plJpa Z[,dji F pa + ^ e p,vprT Z' p d v h (d p Z a - d a Z p ), (3.9) 

where a 2 , a 2 are related to other parameters by gauge invariance as 

a 2 =—— ^ay sin 9]y + ay cos #vp) , at 2 = -^ 0:3 sin Ow + 0:3 cos Ow^j ■ (3.10) 
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The partial decay widths of the CP-odd vector into Z'y, ZZ , W + W , Ay, hZ, and qq, 
respectively are given as follows: 


T'z'(Z'y) = 


OL^rrv 




967r m 


1 - 


1 + xf' 




r z /(fp+ip-) = 


rz'M = 

i>(^) = IS 

rz'te?) = 


3 ./iTfeZ' 
Z'V 1 

487nn^ 

1-xf 


3/2 


«2 


1 - 


1 — 4:Xg‘ 


+ 

3/2 


1 - 4®f£ 


3/2 


+ 4do 


1 - 


(l + 2:4 


(3.11) 


As can be seen clearly from the gauge invariant higher dimensional operators in Eqs. (3.2) 
and (3.3) and can be checked from the effective gauge interactions in Eqs. (3.6) and (3.9), we 
note that the unitarity cutoff of A ~ lOTeV, implies that ki, 2 j «i, 2 ) aq, 2 ) £* 1,2 < O(10~ 2 ) and 
K3,«3 < 0 ( 1 ). 

For a phenomenological study of the spin-1 resonance, we assume that the higher di¬ 
mensional operators given in Eqs. (3.4) and (3.5) come with pure imaginary coefficients, i.e. 
02 = 04 = 0 and 02 = 04 = 0. Then, we get 

k 3 = % = 0, = 0, K 2 = —-—- (ki sin 9\y T K\ cos 9\y ), (3.12) 

rriz 


for CP-even interactions, and, similarly, 

0:3 = 0:3 = 0, 0:2 = 0, 0L2 


—— (cci sin^vv + oti cos Ow), 
m z 


(3.13) 


for CP-odd interactions. There are two free parameters for SM gauge boson couplings in each 
case, K±, k\ and 04 , dq, respectively. In this case, there are no Ay or hZ decay modes of the Z' 
gauge boson while Z^,ZZ and W + W~ decay modes exist. Therefore, the gauge invariance 
of the higher dimensional operators is crucial in correlating between different decay channels 
of the spin-1 resonance. Turning on small couplings to Higgs, we can maintain the diboson 
resonances as hinted by ATLAS and at the same time have a potential to discover or constrain 
the models with spin-1 resonance further by the decay mode into h'y or hZ. Henceforth, in 
order to explain the ATLAS diboson excess from W + W~ axi&ZZ decay modes, we focus on 
a simple parameter choice with K 3 = £3 = 0 for the CP-even and 03 = 0:3 = 0 for the CP- 
odd. In this case, the ratio between W + W~ and ZZ branching fractions remains constant, 
independent of the remaining parameters for both cases, i.e., BR ^z^^zz) ^ ~ 1-56. In Fig. 
3, we show branching fractions of the CP-even vector as a function of the diquark coupling 
(gx) for above choice of parameters. In addition, we have fixed k\ = 0.01 (considering the 
unitarity bound) to maximize the branching fraction into W + W~ and ZZ , as their partial 
decays widths are proportional to k\ for k\ = K 3 = £3 = bi = 0. Numerically very similar 
results are obtained for the CP-odd vector. 
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Figure 3: Branching fractions of the CP-even vector as a function of gx (diquark coupling) for a 
given set of parameters, = 0.01 and K\ = K 3 = £3 = bi = 0. Very similar results are obtained for 
the CP-odd vector and the CP-odd tensor cases. 


For our numerical study, we set cr = cr = 1 for the CP-even (—cl = cr = 1 for the CP 
odd) and ignore the kinetic mixing and mass mixing. We further set 6* = 0 for the CP-even 
case (bi = 0 for the CP-odd case), leaving «q, k\. K 3 , £3, A and gx for the CP-even, and 
aq, aq, a. 3 , A3, A and gx for the CP-odd, respectively, as relevant parameters. Dependence 
on K3, £, 3 , «3 and A3 are weak, and we set them to zero as mentioned above to make a(h 7) 
and a(hZ ) vanish. Furthermore, we conservatively take ki = 0 = aq, for which <7(^7) also 
vanishes. Turning on non-zero values of ki and oq always reduces the branching fractions 
of the diboson signal. Finally, after setting A = 10 TeV, we show in Fig. 4 the production 
cross sections of the CP-even (left panel) and the CP-odd (right panel) vector bosons in the 
ZZ + W + W~ final state (red solid curves). As the resonance is produced by pp collision, it 
can also decays to the dijet final state. The dark yellow-shaded area is disfavored by ATLAS 
dijet searches [42] and the black dotted curves represent Tz'/mz' = 0.15, 0.1, and 0.05, 
respectively. 

The single production cross section itself is explicitly dependent on the coupling, gx, 
only. However, the decay width changes depending on the rest of parameters, which affect 
the shape of the dijet cross section. Our CP-even vector model is the same as one in discussed 
in Ref. [49] , and we are able to use results there by simply rescaling couplings and branching 
fractions in our parameter space. The blue (solid, dashed, dotted) curves labelled by 10 fb , 
300 fb -1 , and 3 ab~ x represent the projected 95% C.L. exclusion contours for 14 TeV LHC, 
respectively. Unfortunately, this projected sensitivity is not available for other resonances, 
and it is not straightforward to recast the results from Ref. [49] due to different efficiencies. 

We note that as shown in the right panel of Fig. 4, the allowed parameter space requires 


11 




8 TeV, QCD=m z =2 TeV, A=10 TeV 



*1 


8 TeV, QCD=m z =2 TeV, A=10 TeV 



Figure 4: Production cross sections (in fb) of the CP-even (left panel) and the CP-odd (right panel) 
vector bosons in the ZZ + W + W~ final state. The dark yellow-shaded region is disfavored by ATLAS 
dijet searches and the black dotted curves represent Tz'/mz> = 0.15, 0.1, and 0.05, respectively. The 
light green-shaded region represents the allowed space to fit the ATLAS diboson data. The blue (solid, 
dashed, dotted) curves represent the projected 95% C.L. exclusion contours for 14 TeV LHC with the 
corresponding luminosity. 


cii > 0.02, which is close to the unitarity limit. Finally, any reasonable deviation from the 
current choice of parameters would be easily allowed, as long as the corresponding limits can 
be avoided in the final states with Z 7 , hZ, and / 17 . 

4 Spin-2 resonances 

The spin-2 resonance with mass mo couples to the SM particles as graviton does, that 
is, 

, (4.1) 

where T^ v is the energy-momentum tensor. We set the spin-2 resonance to couple to the 
energy-momentum tensor for each SM particle with an arbitrary coefficient, which is gauge 
invariant under the SM gauge groups. The energy-momentum tensor with CP-even interac¬ 
tions to the SM gauge bosons are 

T^v = c,Fj lX F YX v + c 2 F^ x F wx „ + c 3 G^G x „ , (4.2) 

where ci, c 2 , and c 3 are constant coefficients parametrizing the relevant coupling strengths. 
Here, we assumed that the spin-2 resonance couples dominantly to the transverse modes of 
SM gauge bosons [50] while the terms proportional to the metric in the energy-momentum 
tensor vanish under the traceless condition. For a heavy spin-2 resonance with mg S> rri\yz, 
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Figure 5: Similar to Fig. 4 but for the CP-even tensor. 


the gauge boson mass terms can be ignored, even if the spin -2 resonance couples to the 
longitudinal modes of gauge bosons as well [50]. 

The partial decay widths of the spin-2 resonance with CP-even interactions into 77 , Z 7 , 
ZZ, W + W~, and gg [50] are 


r e (77) = 


r S (ZZ) = 4*§(1 - 3*0 + 6 (*g) 2 ) , 

M Z t) = fefg? (1 - *g ) 3 (1 + l *g + J(xg) 2 ) , 

r s m' n- j = 77T : 0 - 4 *S- (1 - + 6(*g-) 2 ) 

Ir e (s9) = !^#. 


c 77 = ci cos 2 0 \y + C 2 sin 2 
Czz = C2 cos 2 + Cl sin 2 9 \y 
Cz 1 = (C2 - Cl) sin(26»iy) (4.3) 
cww = 2 c 2 
C 99 = C 3 J 


where again xf is defined in Eq. (2.3). One may notice that for ci = C 2 , the decay mode, 
Q —> Z 7 , vanishes. We note that the branching fractions of the spin-2 resonance are of 
the similar form as the ones of scalar resonances discussed in Section 2 because the spin-2 
resonance decays through gauge invariant operators composed of field strength tensors. One 
can also suppress the diphoton rate by choosing C 2 = —ci / tan 2 9wi which forces all relevant 
branching fractions to be the same as those in the scalar case. Our parameter scan results 
are summarized in Fig. 5. We also find that production cross section, gg —> Q —>• gg , in the 
demonstrated parameter space was small and therefore, there is no constraint from the LHC 
dijet resonance search. 

On the other hand, there is no counterpart of the energy-momentum tensor for CP-odd 
interactions, but the Lorentz invariance and the gauge invariance dictate the detailed form 
of the interactions. Following the similar step as in the CP-even vector case, the CP-odd 
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interactions of the spin-2 resonance to the SM gauge bosons 2 are given by 


where 


rQ _ F p u 

^CP-odd — ’ 


(4.4) 


= «i e pXpa d x Z v F prj + 01 e pXpa d x Z v (d p Z a - d°Z p ) 

+ {a 2 e pXpa d x W~(d p W a+ - d a w p+ ) + ia 2 d x w-(d p w+ - d x w+) + c.c.) 

d x d u hF pX + ^ d A cU ( 5 m z a - d x z p ). (4.5) 

Here a 2 and a 2 are related to other parameters through gauge invariance as 

a 2 = —— ^ai sin 6 w + Si cos Ow'j i a 2 = —^03 sin 6 w + 63 cos dw') ■ (4.6) 

The operators in T pv are induced from higher dimensional gauge-invariant operators such as 
[D X D U H]^F^ X H, [D X D U H]^F^ X H, [D X D U H]^F^ X H 1 and [D X D V H^F™ X H after electroweak 
symmetry breaking.’’ Therefore, the resulting effective CP-odd interactions of the spin-2 
resonance are of strong similarity to those of the spin-1 resonances as discussed in Section 3. 
Hence, the spin-2 resonance can decay into a pair of electroweak gauge bosons or Higgs bosons 
via symmetry breaking terms in T pu . We note that as in the CP-odd vector case, the unitarity 
cutoff of A ~ lOTeV implies that a± }2 , ci\, 2 < 0( 10~ 2 ) and < 23 , 03 < 0(1). 

We note, however, that diquark CP-odd operator, q'y 5 ('y fl d l/ + 7 v d p )q + h.c., is a total 
derivative, while a nontrivial diquark operator, iq^ d v F'y 1 ' d p )q+\i.c. is CP-even [53]. Since 
the CP-odd diquark operator can be written as (d v G p , v )q'y b 'y p, q by integration by parts, the 
diquark production of the on-shell CP-odd spin -2 resonance is suppressed due to d u G pv = 0. 
Instead, the CP-odd spin-2 resonance can be produced by vector boson fusion. In this case, 
there are two forward jets accompanying the resonance, so we cannot explain the ATLAS 
diboson excess by the CP-odd spin-2 resonance. For this reason, we do not consider it 
any longer. We also remark that the operators composed of field strength tensors only, for 
example, G pu Tv(F pX F x u ) with F /w being F pu , F™, or Gjju , vanish because G Tr ( F !lX F ^) — 
— \Gjl Tr(F a pF al3 ) = 0 due to the traceless condition, i.e., G/l = 0. Therefore, those gauge 
invariant operators do not contribute to the process with on-shell CP-odd tensor so that the 
gluon fusion production of the CP-odd spin-2 resonance is suppressed. For a future reference 
on the phenomenological study of the CP-odd spin -2 resonance, we summarize the partial 

J The ZZ coupling to the CP-odd tensor field was considered in Ref. [51, 52] without gauge invariance 
imposed. 

j We note that one of higher dimensional operators among [D V D X HyF X X H and [D x D u Hy FFH is redun¬ 
dant because [D X D V H - D v D x H] f F Y x H ~ \H\ 2 F YX V F Y X , which contributes to the gauge invariant operators 
of F yx uFlx that becomes a vanishing gauge interaction after electroweak symmetry breaking (EWSB). 
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decay rates with CP-odd interactions into Z7, ZZ, W + W —, /17, hZ and qq as listed below: 
r 5 (^7) = S (1 - 4) 3 (34 + 3*g + 3(*f )- 1 ) , 

r ^(^) = \A " ^ ( 3 " “ 32( ^ )2 ) ’ 

r g (W+W~) = 1920 ^ a2 (1 - 4x^) 3/2 [3a| (1 - 4x% ) + 4a 2 (3 + 8*® )] , 

^(^7) — 1280 ttA 4 (■*■ — ) ’ , (4.7) 


r 0{hZ) = 3& 


1 - x h 


- 2 ( 1 + x z) x z + 04 ') 


„G\2 


5/2 


s(i - x z) + 2 ( - 2 + 5x£ + 04) 2 )zf - (l + 4xf) (xg ) 2 + 2(xf ) 3 


1280ttA 4 


where again xf is defined in Eq. (2.3). 


5 Kinematic Correlations in the Diboson Final State 


In this section, we discuss ways of discriminating potential scenarios to give rise to diboson 
resonances. Since we have observed that various bosonic particles with different spins and CP 
states can accommodate the excesses reported by the ATLAS collaboration with a suitable 
choice of parameters, it is of paramount importance to pin down the underlying physics 
once those excesses are confirmed experimentally. Of potentially useful variables, we employ 
several angular correlations between the decay products of the resonance of interest. We first 
suppose that a resonance R decays into two vector bosons Vi and F 2 which subsequently 
decay into two visible particles u, and x* (* = 1 , 2 ): 


pp —> R —> Vi(—> Ui + V\) + Ih 2 (—> u 2 + V 2 ), 


(5.1) 


and denote Pi as the three momentum of Vj and Hi(vi ) as those of Ui(vi). 

With these notations, we enumerate the angular variables to be used here as follows: 


Pi • (hi x h 2 ) 
Pi ■ (hi x h 2 ) | 


cos 


(hi • h 2 ) with hi 


Ui X Vj 
Ui X Vi \’ 


+ 1 


p } ' X hsc) cos - 1 (hi • h sc ) with h sc 
Pi ■ (hi x h sc )\ 


z x Pi 
z x Pi 


cos 9* 


Pi-z 

\Pi\ 


COS@l = 


P2 • Ul 
\P 2 \\UI 


(5.2) 

(5.3) 

(5.4) 

(5.5) 


For the first three variables, all the momenta are measured in the rest frame of resonance R, 
while for the last one, all the momenta are measured in the rest frame of vector boson V\. 
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A.U A.U 


Scenario 

Parameter choice 

R production 

0 + 

si = 0.4, S 2 = —si/tan 2 S 3 = 1, A = 10 TeV 

99^ R 

O' 

ai = 0.6, 02 = —a\/ tan 2 Q\y, 03 = 1, A = 20 TeV 

99^ R 

1 + 

k\ = 0.008, gx = 0.02, cl = cr = 1, A = 10 TeV 

qq —>• R 

1 ' 

di = 0.01, gx = 0.04, — cl = cr = 1, A = 10 TeV 

qq —>• R 

2 + 

ci = 0.5, C 2 = —ci/ tan 2 0p/, C 3 = 1, A = 5 TeV 

99 -t R 


Table 1: List of scenario choices for a resonance R having a spin and CP-state denoted as J CP . 




* *! 



Figure 6: Unit-normalized distributions in d> (upper left panel), $1 (upper right panel), cos0* (lower 
left panel), and cos6 \ (lower right panel) for the resonance decay into two W gauge bosons. The spin 
and CP state of the resonance of interest is represented by J CP . 


These variables have been used in the context of resonance discrimination in the literature, 
and they show distinctive structures depending on quantum numbers of each resonance (see, 
for example, Ref. [51]). 
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We here study the above-listed observables in the analysis of R —> W + W~, and show 
the distributions in Fig. 6 . The distributions are plotted with parton-level events with a 
10% of Gaussian smearing onto energy of each final state particle for more realistic Monte 
Carlo simulation. Again, events were generated by MadGraph_aMC@NLO [54] together with 
the default set of parton distributions NNPDF23 [55] at the center of mass energy of 13 TeV. 
Table 1 summarizes our parameter choices for each scenario. All the parameters not listed 
in the table are simply taken to be zero. Note that this choice of parameters is made only 
for the purpose of illustration of different kinematic distributions for each scenario. We find 
that the shape is not strongly dependent on parameters. The mass and the total decay width 
are fixed to be 2 TeV and 0.1 TeV, correspondingly. We remark that the spin-0 and CP-even 
spin -2 resonances are produced via gluon fusion while the spin -1 resonances are produced via 
quark annihilation. The observables of <b, <Fi, cos#*, and cos#i are exhibited in the upper 
left panel, the upper right panel, the lower left panel, and the lower right panel, respectively. 
Different spin and CP states are symbolized by J CP , and they are histogramed as follows: 
CP-even scalar by the blue dashed, CP-odd scalar by the red dashed, CP-even vector by the 
green solid, CP-odd vector by the orange solid, and CP-even tensor by the black dot-dashed. 
In particular, the theory prediction for cos 6 * distributions is readily derived as follows: 

1, for gg ^ 0 + , 0 —>• W + W~ 

1 +cos 2 #*, for qq -> 1 +, 1" -> W+W~ . (5.6) 

1 + 6 cos 2 6 * + cos 4 6 * , for gg -> 2 + -> W+W~ 

which can be directly compared with experimental data. 

First of all, we observe that the angular distributions with the Gaussian smearing are very 
similar to those without any smearing, from which we expect that the angular distributions 
are insensitive to detector effects such as jet energy resolution. Moving onto Fig. 6 , we clearly 
see that these observables are useful enough to distinguish potential scenarios associated with 
diboson resonances. For example, the CP-even vector resonance (green solid histograms) 
shows distinctive behaviors in all four observables. Furthermore, the unique features according 
to different spin and CP states in those variables can be used for cross-checks. Note that 
one single distribution can not discriminate different scenarios, and thus it is important to 
consider all possible kinematic correlations. Finally, we remark that similar analyses can 
be straightforwardly applicable to other diboson resonances such as R —> h'y, R —> hZ 
and R —> Zq so that more information can be extracted to confirm the underlying physics 
governing the observed phenomena. 

6 Summary 

Recently, the ATLAS collaboration has reported some excesses in searches for diboson res¬ 
onances using jet-substructure techniques. The excesses show up in the invariant mass of 
W + W~, W P Z and ZZ at around 2 TeV. It has been discussed in literature that about 20% 


da 

- rsj < 

d cos 6* 
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of the events in at least one signal region belong to all three categories, which indicates that 
these “resonances” may be explained by one single particle rather than two. 

In this paper, we have explored a possible new physics interpretation of the ATLAS 
diboson excess in an effective field theory approach, which covers a rather large class of 
models in a reasonably model independent manner. We considered the effective operators for 
scalar (s = 0), vector (s = 1), and tensor (s = 2) resonances with different CP properties. It is 
shown that each scenario may explain the ATLAS diboson excess without contradicting other 
constraints, except the CP-odd spin-2 resonance whose diquark or gluon fusion production 
is suppressed. The CP-odd vector case might have some tension with the unitary bound. 
Symmetries of each scenario predict signals in other final states such as Zj and 77 in the 
cases of scalar and CP-even tensor resonances; Z 7 and hZ, /17 at a smaller rate in the cases 
of vector resonances. Especially, the dijet, tt, Z 7, hZ, and /17 resonance searches at the LHC 
run II may confirm or constrain these scenarios. 

With limited statistics, all these scenarios provide a relatively good fit to the data. How¬ 
ever, a further accumulation of data might reveal the real identity of the resonance. We 
showed a few examples of kinematic distributions, which are sensitive to the CP property 
and spin of the resonance. We strongly encourage experimental collaborations to look at 
these kinematic correlations. 
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A Decay widths 

In this appendix, we summarize the useful formulas for the decay rates for scalar and tensor 
resonances. 

A.l CP-even scalar 

For interaction Lagrangian C = — cy, y 2 ^ Fy'Fy 2 /t „. the decay width of 4> to V\ Vj is given as 

(A.i) 

87 r y A- J 171 ^ m$ 
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where F(x\,x 2 ) is defined as 


T{x\,x 2 ) = (l — (xi + x 2 ) 2 ) 1 ^ 2 (1 - (xi - x 2 ) 2 ) 1 ^ 2 (l + xf + x\ — 2(.x 2 + x|) + 4 x 2 X2) . 

(A.2) 


sy is symmetric factor, which is 1 for V\ / V 2 and 2 for V\ = V 2 , respectively. 


A.2 CP-odd scalar 

For interaction Lagrangian C = — cy y 2 ^ Fy 2 /t „, the decay width of A to F 1 V 2 is given as 


T(A^ FiF 2 ) = 


s v\cv 1 v 2 1 
2vr 


m 


mi m 2 




t/(xi,x 2 ) = (1 - {xi +x 2 ) 2 ) 3/2 (1 - (xi - x 2 ) 2 ) 


m A m A 

x2\3/2 


(A.3) 

(A.4) 


where sy is symmetric factor, which is 1 for V\^V 2 and 2 for V\ = V 2 , respectively. 


A.3 CP-even tensor 

For spin-2 tensor with mass me, the interaction Lagrangian is C = — c v 1 V 2 ^x 2 Fv 1 \Fv 2 Xl ' and 
the decay width of Q^ u —»• V\V 2 is given as 


r (g^ -»V1V2) 


sv\c Vl y 2 \ 2 m 3 G 

1607rA 2 


m 1 
m h ’ 


m 2 


) 5 

m h 


(A.5) 


where the convenient dimensionless function, T-L(x,y), for some interesting cases are 


T-L(x, x) = \/l — 4.x 2 (1 — 3x 2 + 6x 4 ), (A.6) 

W(x, 0) = (1 - .x 2 ) 3 (l + \x 2 + |x 4 ) (A.7) 


and the symmetric factor sy is 1 for V\^V 2 and 2 for V± = V 2 , respectively. 
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